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ABSTRACT 

Context. The physics of galaxy clusters has proven to be influenced by several processes connected with their galactic component 
which pollutes the intracluster medium (ICM) with metals, stars and dust. However, it is not clear whether the presence of diffuse 
dust can play a role in clusters physics since a characterisation of the infrared (IR) properties of galaxy clusters is very challenging 
. and yet to be completely achieved. 

Aims. In our study we focus on the recent work of Giard et al. (2008) who performed a stacking analysis of the IRAS data in the 
direction of several thousands of galaxy clusters, providing a statistical characterisation of their IR luminosity and redshift evolution. 
We model the IR properties of the galactic population of the SDSS-maxBCG clusters (0. 1 < z < 0.3) in order to check if it accounts 
O . for the entire observed signal and to constrain the possible presence of other components, like dust in the ICM. 

Methods. Starting from the optical properties of the galaxies of the SDSS-maxBCG clusters, we estimate their emission in the 60 and 
100 /jm IRAS bands making use of modeled spectral energy distributions of different spectral types (E/S0, Sa, Sb, Sc and starburst). 
We also consider the evolution of the galactic population/luminosity with redshift. 
Results. The total galactic emission, which is dominated by the contribution of star-forming late-type galaxies, is consistent with the 
observed signal. In fact, our galactic emission model slightly overestimates the observed fluxes, with the excess being concentrated in 
low-redshift clusters (z < 0.17). 

Conclusions. Our results indicate that, if present, the IR emission from intracluster dust must be very small compared to the one 
associated to the galaxy members. This translates into an upper limit on the dust-to-gas mass ratio in the ICM of Z<j < 5 x 10~ 5 . 
The excess in luminosity obtained at low redshift constitutes an indication that the cluster environment is driving a process of star- 
formation quenching in its galaxy members. 
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, 1. Introduction AGNs feedback and other non-gravitation al processes, such as 



supernovae (SNe) p owere d galactic winds (iKapferer et al. ll2006 ; 
Sijacki & Spri ngeill2006t jSchindler & Diaferiol|2008l) . preheat- 



Clusters of galaxies form in correspondence of the peaks of the 

primordial matter density field as a result of the gravitational col- ing of the ICM (see e.g. iFang & Haimanll2008l) . together with 

lapse of both dark matter and baryons. In the framework of the gravitational ones (i.e. galaxy-galaxy interactions, ram-pressure 

standard ACDM cosmological model and the hierarchical clus- stripping, galaxy mergers) induce energy and matter exchanges 



X 

^ , tering of the large scale structure formation they constitute both between the galactic medium and the ICM. To date these com 

the most recent and the largest virialised objects of the Universe, plex physical processes and their impact on the statistical clus- 

Nowadays, it is clear that gravitation is not the only pro- ter PjTOtie^thus on our understanding of structure forma- 

cess that influences the physics of the intracluster medium (ICM, tlon (|Voitl|2005|) and use of the cluster population in cosmologi- 



hereafter). The electrons of the ionized plasma emit via free- cal studies (|Mantz et al.|2008|;|Vikhlinin et al.|2009|:|Mantz et al.| 



free interaction with the protons, making clusters bright X-ray EQOjh are under serotinous observational and theoretical inves- 

sources and allowing the gas to cool efficiently, particularly in ligations (see, e.g., |Arnaud|2005t |Borgani|20084 and references 

the central regions. During the last decade, the observations of therein). 

the XMM-Newton and Chandra X-ray satellites highlighted the The presence of heavy elements in the ICM is the most di- 
presence of several interaction mechanisms between the galactic rect evidence and consequence of the ejection of galactic ma- 
component and the ICM, showing that the evolution of the two terial. Within clusters their abundance has been widely mea- 
is intimately tied. For instance, the accretion of cold gas onto the sured making use of X-ray observations (see the reviews by 
brightest cluster galaxies (BCGs) at the cluster center is strongly ISarazi rnau d 2005^ IWerner et al.ll2008h with a typical 
suspected to fuel the super-massive blackholes they host. This abundance of 0.3 Z . Stars and SNe constitute the most effi- 
process is able to trigger star-form ation within the BCGs cient way to produce and disperse metals. Heavy elements in 
dO'Dea et al.l l2008; Pipi no et alj 12009) and power episodic vi- the ICM originate from different processes: early enrichment 
olent outbursts of their central active galactic nuclei (AGNs), (lAguirre & Schavel 120071) . continuous injection from galaxy 
whose energy injection into the ICM prevents the overcool - members and in situ product ion by intra-cluster stars (sources 
ing of the gas (Fab ian et al.ll200oT: iMcNamara & Nul sen 2007). of the intra-cluster light, see lArnaboldil2004l: iRrick & Bernstein 
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20071: iMurante et all l2007t IConrov et alj 120071: iDolag et alj 



2009). The aforementioned processes do not discriminate be- 



tween the natures of the ejected galactic material, therefore these 
enrichments are unambiguously linked also to the ejection of 
neutral gas and dust in the ICM. 

Observations indicate that in galactic environments dust is 
a minor component , with dust-to-gas ratio Mdust/^gas ~ 0.01 
(Mat his etal.lll977|) and it could be as low as Md us t/M gas = 
10~ 5 - 10~ 4 in the ICM dPopescu et alj l2000t lAguirre et al.1 
1200 ll) . Nonetheless, dust particles have lifetimes long enough to 
be heated by collisions with the hot electrons and re-emit at the 
infrared (IR) wavelengths. In this way they can constitute an ad- 
ditional cooling agent of the gas dMontier & Giardll2004l) which 
might play an important role in ICM physics, as seen with the 
implementation of d ust cooling in hydrodynamical simulations 
(Ida Silva et al.ll2009h . 

However obtaining observational constraints on the possible 
IR signal coming from intracluster dust is a very challenging 
issue, since the average sky fluctuations caused by background 
galaxies and galactic cirrus clouds are comparable or even higher 
than the overall flux coming from a single cluster, which is any- 
way expected to be dominated by the dust emission from star- 
forming galaxies. In fact, nowadays the only claimed (and still 
controversial) detection of this IR emiss i on comes fro m the stud- 
ies on the Coma cluster o f lStickel et al1(ll998ll2002l) . who mea- 
sured a diffuse dust mass of Md us t ~ 10 7 M(n, while othe r at- 
tempts resulted in non detections (see, e.g., Bai et al. 2007). 

If studies on single objects appear very problematic, the low 
signal-to-noise problem can be overcome by adopting a statis- 
tical approach, taking advantage of the high n umber of clus- 
ters detected mainly with optical surveys (see iBivian ol 12001 
for a review). Th i s conc ept was first used at IR wavelengths by 
Mont ier & Giardl (|2005) who performed a stacking analysis of 
the Infra Red Astronomical Satellite (IRAS) survey in the di- 
rection of more than 1 1 000 known galaxy clusters and groups 
and obtained a statistical detection of the overall clu sters signal 
at 12, 25, 60 and 100 p.m. Starting from this result. IGia rd et alj 
(2008) characterised the IR luminosity evolution of the stacked 
sample, analysing also its correlation with their X-ray luminosi- 
ties. These works constitute a first important result on the global 
IR properties of galaxy clusters. However, since the clusters IR 
emission could be dominated by star-forming galaxies, it be- 
comes crucial to disentangle the galactic signal from a possible 
diffuse component in order to quantify its implications on clus- 
ter physics. On the other side, understanding the IR properties 
of the cluster galaxies could provide hints on the environmental 
effects over their star-formation, which are nowadays objects of 
studies a t different wavelengths, fro m the radio to the ultraviolet 
(UV, see iBoselli & Gavaz zi 2006a, for a review). 

In this work we attempt to reconstruct the observed stacked 
IR emission measured in the direction of galaxy clusters in 
the IRAS whole sky survey. We model in details the contri- 
bution of cluster galaxies to the total IR emission in order to 
understand whether this contributi on is sufficient to explai n the 
fluxes and lumino sities derived bv lMontier & Giardl d2005l) and 
iGiard et aD d2008l) or if, on the contrary, there is an indication 
of a non-galactic component, possibly associated to intraclus- 
ter dust. Since for our modelisation we make use of the known 
spectral properties of galaxies as observed mainly in the field, 
the comparison between our results and observed data can be 
useful also to highlight possible differences between cluster and 
field galaxies. 

This paper is organised as follows. In the next Section we 
briefly present the statistical detection of the clusters IR emis- 



sion of iMontier & Giardl d2005l) and IGiard etail d2008l) and we 
focus in particular on the subsample of the SDSS-maxBCG cat- 
alogue. In Section [3] we describe the details of our modelisation 
for the IR luminosity of cluster galaxies and discuss the possible 
presence of other galactic components in Section|4] In Section|5] 
we apply our model to estimate the IRAS fluxes, then our results 
are discussed in Section [6] Finally, we present our conclusions 
in Section|7] 

Throughout this paper we assume a flat ACDM cosmological 
model with (Cl A ,£l m ,h) = (0.7,0.3,0.7). 

2. The statistical IR emission of clusters 

Working on a list of 1 1 507 gr oups and clusters se l ected from 
14 publicly available catalogues^ Montier & Giardl d2005l) per- 
formed a clear statistical detection of the IR flux in the direc- 
tion of galaxy clusters by stacking their corresponding fields 
(within 10' from the cluster center) in the IRAS all sky sur- 
vey. After dealing carefully with the point sources contamina- 
tion, foreground/background subtraction and other various sys- 
tematic effects, they measured stacked fluxes at 60 and 100 /jm 
with signal-to-noise ratio of 57 and 43 respectively. Over the 
four IRAS wavelengths this emission proves to be consistent 

with the spectral signature of galac tic IR emission. 

More recently, on th e basis of IMontier & Giardl (120051) re- 
sults, |&ardetal] d2008) performed a statistical analysis in red- 
shift, presenting for the first time the evolution of the IR lumi- 
nosity of galaxy clusters. In additi on to IMontier & Giardl (f2005) 
original list of 1 1 507 clusters, IGiard et all d2008l) backed-up 
their analysis on two standalone catalogues - i.e. the "Northern 



Sky Optical Cluster Surve y " (NSC. iGal et alJ 120031) . and the 
SDSS-maxBCG catalogue dKoester et al.l 120071) - for which a 
richness information was available. They extended their IRAS 
stack ing analysis to the Rosat All Sky Survey (RASS, [Voges 
1 19921) in order to compare the evolution in redshifts of the IR 
and X-ray luminosity. They showed that the stacked IR lumi- 
nosities are on average 20 times higher than the X-ray lumi- 
nosities. They also found that the IR lumi nosity is evolving 
rapidly as (1 +z) 5 in the 0.1 < z < 1 interval. IGiard et al.l(l2008l) 
also made use of the richness information (i.e. correlated to the 
halo occupation number), contained in the SDSS-maxBCG and 
NSC catalogues, to constrain the dependence of the IR luminosi- 
ties with cluster richness. They derived a correlation following 

In order to understand the aforementioned results at 60 and 
100 /jm, in this work we focus on the modelisation of the IR lu- 
minosity and fluxes of galaxy clusters. The I R spectrum of galax- 
ies, which i s dominated by dust emission (Lag ache et al.l l2005: 
ISoifer et al.l [2008 ). is expected to make a major co ntribution to 
the IR emission detected in the direction of clusters (Giar d et aD 
I2008I Therefore, a careful modeling of the IR emission of mem- 
ber galaxies is needed in order to understand the global proper- 
ties of clusters in the IR, and in particular to disentangle between 
the relative contribution of galaxies and intrac luster dust. Given 
the inh omo geneity of the main dataset used bv lMontier & Giardl 
(120051) and lbiard etail d2008l) (i.e. the 1 1 507 groups and clus- 
ters) in terms of cluster size and detection method, it becomes 



The list was build making use of the SIMBAD database, oper- 
ated at CDS, Strasbourg, France (http://simbad.u-strasbg.fr/simbad/ 1 
and NASA/IPAC Extragalactic Database (NED), operated by the 
JPL/Caltech (http://nedwww.ipac.caltech.edu I, among which the 
ABELL catalogue (Abelletal. 1989) and the "Northern Sky Optical 
Cluster Survey" (NSC. lGal et alj|2003l) . 
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very difficult to define a global selection function of such sam- 
ple, that would be necessary in order to have a suitable starting 
base for our work. In order to overcome these difficulties we de- 
cided to base our modelisation on a restricted and well defined 
cluster sample by considering a single catalogue. For the purpose 
of our calculations we chose the SDSS-maxBCG catalogue. 



2.1. The SDSS-maxBCG catalogue 

Since we want to reconstruct the total IR luminosity and flux 
due to cluster galaxies, we require an information on the halo 
occupation number (i.e. the number of galaxies within the clus- 
ter potential well). In this view the SDSS-maxBCG catalogue 
( Koes ter et al.l2 007) is the most fitted to our needs. Indeed, since 
the SDSS-maxBCG catalogue is created by identifying overden- 
sities in the galaxy distribution, it contains information about 
the richness of the cluster itself that we can use directly without 
having to deal with the uncertainties connected with the N — M 
and Lx — M relations, that would be necessary, for instance, in 
case of X-ray selected clusters. Moreover, the number of objects 
contained in the SDSS-maxBCG catalogue is the highest among 
all the others (i.e. 13 807 groups and clusters), thus providing a 
good statistics which is needed to a ssess the robustness of our 
results. Finally, iKoester et alJ d2007h provide also the total lu- 
minosity in the r and /-band of the identified cluster members 
which also have a well defined spectral type (E/SO galaxies). We 
will make use of this information to obtain an estimate of their 
flux in the IRAS bands (see Sect.0. 

The SDSS-maxBCG catalogue has been obtained by 
analysing the clustering properties of more than 500 000 SDSS 
galaxies in an area of ^7 500 deg 2 in the redshift range 0.1 < 
z < 0.3. For each cluster the catalogue contains a photometric 
redshift and an indication of the cluster richness, A^j 00 . This 
quantity corresponds to the number of early-type galaxy mem- 
bers at a distance lower than /?200 from the central BCG The 
adopted definition of R200 is the radius at which the deprojected 
galactic density is 200 £2 m ~' times the average galactic density 
on the large scale: in the approximation of galaxies following 
the matter distribution, this is equivalent to the usual cosmologi- 
cal definition of R200 as the radius enclosing a density 200 times 
la rger than the critical density of the Universe (see the discussion 
in lHansen et al . 2005, for details). 

It is important to note that A^J does not represent the to- 
tal cluster richness because it does not include spiral galaxies: 
in fact, the cluster members candidates hav e been selected for 
having colors matching the E/SO ridgeline dBower et alJll992l) 
and M r < — 16. We will explain in the detail our definition of the 
cluster richness in Section [3~2l 

The SDSS-maxBCG catalogue includes 13 823 groups and 
clusters with N^°> 10 and it is widely dominated by the pres- 
ence of groups and small clusters (N™°~ 10 — 15). As a refer- 
ence, when computing the ir corresponding M200 with the rela- 
tion of Rvkoff et al. ( 2008), the median of the mass distribution 
is 4.5 x 10 13 h~ l M P) . I n this work, we adopted the trimming done 
bv lGiard etal] d2008l) on the SDSS-maxBCG catalogue of ob- 
jects located in regions of the sky not covered by the IRAS data, 
with noisy images, or whose fields contain strong IR sources. 
This selection process ended up in a sample of 7 476 clusters, 
corresponding to a total of 121 318 E/SO galaxy members, that 
will be the object of our analysis. 
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Fig. 1. SEDs as a function of wavelength for the 5 GRASIL 
dSilva et al.l 119981) templates used in this work. The red line 
shows our reference model for E/SO galaxies, the three blue lines 
represent normal spiral galaxies Sa, Sb, Sc (light blue, average, 
dark blue, respectively) and the green line represents our star- 
burst galaxy model, corresponding to M82. All of the templates 
have been arbitrarily normalized to the same integrated luminos- 
ity of 10 30 erg/s in the band 8-1000 /urn (i.e. IR luminosity, Lir). 
The four shaded regions identify the optical r and /-band and the 
two IRAS bands of 60 and 100 /jm, from left to right. 

3. Modeling galaxy luminosities 

The main source of the IR radiation of galaxy clusters is ex- 
pected to be the thermal emission due to the dust reprocessing 
the UV photons emitted by stars inside cluster galaxies. Since we 
are basing our work on galaxies observed in the optical band, we 
need to construct a valid model that is able to connect their galac- 
tic r band emission to the corresponding dust emission in the 
IR. In particular, we will focus on the 60 and 100 /jm IRAS 
bands fluxes, in order to compare our re sults with the most rele- 
vant measurements of iGiard et al.l d2008l) . The characteristics of 
the galactic IR emission depend on the amount of dust present in 
the galaxies and are strictly connected with their star-formation 
rate (SFR) and, therefore, their spectral type. 

Spectral energy distribution (SEP ) templates ob tained with 
the spectral synthesis code GRASIL dSilva et alJ [l998) were used 
in order to represent the typical emission of galaxies in the local 
universe. GRASIL SEDs extend from UV to radio wavelengths, 
including dust reprocessing and nebular line emission. We used 
four templates for normal galaxies, with Salpeter's initial mass 
function and age 13Gyr, namely an elliptical galaxy and three 
different spiral galaxies (Sa, Sb and Sc). We also included the 
GRASIL model fit to multi-wavelength observations of M82 as 
a semi-empirical templ ate to represen t the typical SED of lo- 
cal starburst galaxies dSilva et al] [l998). For spiral galaxies, the 
SEDs correspond to a weighted average over the different lines 
of sight, from face-on to edge-on, in order to statistically ac- 
count for the mean spatial orientations of cluster galaxies. The 
five templates are displayed in Fig.Q] 

The SED models described above were used to derive the ex- 
pected IR luminosities from observed fluxes and corresponding 
luminosities in the r band. For a given template, the luminosity 
in the r band is given by 

Lr= f L(X)T r (X)dX , (1) 
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Table 1. Reference spectral types and galactic population mod- 
els used in this work. 



Spectral type 


Population (%) 


u — r 






no ev. 


ref. 


cons. 






-Early-type 












BCG 


4.3 


4.0 


4.2 




0.12 


E/SO 


65.4 


61.5 


63.3 




0.12 


(Tot. early-type) 


(69.7) 


(65.5) 


(67.5) 






-Late-type 












Sa 


4.4 


5.1 


6.5 


>2.2 


2.36 


Sb 


9.6 


10.7 


11.7 


1.8-2.2 


13.82 


Sc 


15.8 


18.2 


14.1 


1.0-1.8 


33.80 


Starburst 


0.5 


0.5 


0.2 


<1.0 


73.16 


(Tot. late-type) 


(30.3) 


(34.5) 


(32.5) 







The fraction (in percentile) of each type of galaxies in the three pop- 
ulation models is shown: without redshift evolution, in our reference 
model and for the conservative scenario (second to fourth column, re- 
spectively). The fifth column shows the corresponding intrinsic u — r 
color intervals used for the modelisation of the spirals (early-type galax- 
ies are taken directly from the SDSS-maxBCG catalogue). The last col- 
umn shows the values of the ratio Lir/L,- for each template. 

where L(X) is the luminosity per unit wavelength and T r (X) is 
the SDSS r filter transmission. For the purpose of this wo rk, we 
define the IR luminosity following lLe Floc'h et al. (2001: 



/.lOOOj/m 

L IR = / L{X)dX. 



(2) 



Therefore, for each template SED, the scaling ratio can be de- 
fined as follows 



K lR,r = -j— 

Lr 



(3) 



Table Q] summarizes the R\^ r values for the different tem- 
plates used in this paper. For a given L r , the corresponding Lfr 
can be different by almost three orders of magnitude depend- 
ing on the spectral type of the galaxy. The modelisation of the 
luminosity of early-type and late-type galaxy populations is ex- 
plained in the details in the next sections. 

3.1. Early-type galaxies 

Early-type galaxies are characterized by an old stellar popula- 
tion and a star-formation history which is essentially compatible 
with passive evolution: for this reason they are not expected to 
dominate the emission in the IRAS bands. In fact, for a fixed 
Lir, the amount of energy emitted at wavelengths X <; 40 
is one order of magnitude lower with respect to normal spirals 
(see FigJTJ. However, it is known that in dense environments 
elliptical galaxies largely dominate the galactic popul ation, be- 
ing about 4 times mo re frequent than spirals (see, e.g.,KSi 
1980; Dressier et al.1ll997l) . Moreover red galaxies are usually 
more massive and more luminous in the r band than blue ones. 
For these reasons their total contribution on the clusters IR signal 
may be non-negligible. 

As mentioned in Section lXTl Koest eret al.l (120071) identified 
the number A^ 00 of E/S0 galaxies (M r < -16) inside fl 2 oo of 
each of the 7 476 clusters of our sample. For each cluster, the au- 
thors provide the luminosity, ^-corrected at z=0.25, in the r band, 
of the BCG and of the other E/S0 members as a whole. We cor- 
rected these luminosities into rest-frame luminosities by using 



the LRG template of KCORRECT dBlanton & Roweisll2007l) : we 
refer to L? CG and L™ emb for the luminosity of the BCG and of the 
cluster as a whole, respectively, after this correction. Given the 
uniform properties of early-type galaxies, and the high number 
of such objects in our sample, their typical SED should be well 
represented by the E/S0 template described above. In fact, al- 
though the IR signal of every single galaxy is lost in the stacking 
process, the E/S0 template still represents the average behavior 
of the early-type population of galaxies. 

For every cluster we assign a r-band luminosity to all of its 
early-type members contained in the SDSS-maxBCG catalogue. 
For the BCG we use the value of L® CG , while for the other galax- 
ies we use the average luminosity of the non-BCG E/S0 galaxies, 
determined as 

r memb j BCG 



£ avg = 



V™ 



N; 



R2Q0 



gal 



1 



(4) 



From eq. [3] we translated these r-band luminosities into IR lu- 
minosities: Lir = RjR,r x L r , subsequently used to normalise the 
SED and further on to derive the fluxes in the 60 and 100 /jm 
IRAS bands (see Sect.PJ). 

3.2. Contribution of late-type galaxies 

Although the majority of optically bright galaxies in clusters 
environment are elliptical, and despite spiral galaxies in dense 
environments tend to be quickly stripped of their gas and have 
their star-formation quenched, they are still expected to provide 
a dominant contribution to the IR emission, due to their higher 
SFRs. Therefore, the contribution of late-type galaxies to the IR 
emission is crucial for our purposes and should be carefully es- 
timated. Given the fact that late-type galaxies are not considered 
in the SDSS-maxBCG catalogue, their contribution was com- 
puted based on available modelings of galaxy populations, in 
particular the distribution of galaxies as a function of spectral 
type and luminosity. Since our sample goes from galaxy groups 
to rich clusters, we have included the environment dependence 
of these variables. 

We based our model on the known prope r ties of galaxies in 
the local universe as given by iBalogh et al.l (120041) who anal- 
ysed the bimodal distribution of galaxies in a local (z < 0.08) 
sample of SDSS galaxies (DR1). These authors performed a de- 
tailed study on the color and luminosity distribution of galax- 
ies as a function of their environment for both early and late- 
type galaxies. For each galaxy they define a density estimator £5 
which represents the local projected galaxy density (M r < —20) 
in Mpc and they use it to classify the different environmen- 
tal regimes. Beside mentioned before (see Sect. 12. lb . 
the SDSS-maxBCG dataset contains also information about the 
number of galaxies N ga \ projected within a distance lower than 1 
h~ l Mpc from the BCG, with the magnitude limit of M r < -16. 
Therefore for each cluster we can assume a unique value of £5 
for all of its members by establishing a relation with iV ga i and, fi- 
nally, with the properties of the population of late-type galaxies 
present. 

We proceed as follows. 

i) First of all, we integrate the do uble gaussian distri butions of 
red and blue galaxies (Fig. 1 of lBalogh et al . 2004) obtaining 
a ratio of red and blue galaxies for the 5 values of £5 con- 
sidered: for our work, we only need the three density bins 
corresponding to dense environments (i.e. groups and clus- 
ters, £5 > 0.5). In these bins the resulting spiral fractions are 
/spi =0.597, 0.478 and 0.286, in order of increasing density. 
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11) 



iii) 



In order to use these values for our calculations we need, 
at first, to correct th e value of N ga \. The galaxy samples of 
iKoester et al.1 (120071) . in fact, has a magnitude limit M r < 
— 16 while the definition of £5 includes only galaxies with 
M r < —20: for this reas on we use the red gal axies luminosity 
function (LF) given bv lBaldrvet all (2004) to calculate the 
number of expected galaxy members with M r < —20. We ob- 
tained < N~^° >=0.32jV sa i. A s for lBaloghetal] d2004l) . the 
sample of Bal drv et alj (|2004) is constituted by local SDSS 
galaxies, with the difference that it contains no density dis- 
tinction (e.g . field galaxies). Therefore, by using the LF of 
iBaldrv et al.l (J2004) for this calculation we neglect the effect 
of environmental properties; anyway, since red galaxies are 
mostly present in galaxy clusters, we do not expect this to 
change significantly our results. 

We are now able to obtain a first estimation of the density 
parameter associated to every cluster with the formula 



£5,0 = 



1-Zspi.O rc(lMpc) 2 



(5) 



iv) 



v) 



where / S pi,o is me fraction of spiral galaxy (with M r < —20) 
that for this calculation we fix arbitrarily to the initial value 
to/spi,o = 0.3. 

Using the value of £5 we calculate a new estimate /spi,i of 
the spiral fr action by interpola ting between the values of / sp j 
obtained by Balogh et al. (2004); this allows the calculation 
of a new value £5 \ of the galaxy density with equation[5] 
We repeat this procedure until E5 , and / sp i.,- converge to def- 
inite values and check that the result does not depend on the 
initial choice of / sp i,o, thus obtaining the expected number of 
spiral galaxies with M r < —20 



< AT 20 >=< A' 



-20 

gal 



> X- 



spi 



N: 



R2QQ 



gal 



spi 



N. 



(6) 



gal 



where the factor Nj? 200 ' 



y , /iVgai accounts for the ratio between 

the galaxies inside R200 and \hr x Mpc. 
vi) Finally, we add the expected number of spiral galaxies in the 
magnitude range —20 < M r < — 1 8 by refer ring again to the 
galaxy distributions of Balo gh et alj ((2004), thus obtaining 
for every cluster the estimated value of the number of spiral 
members with M r < — 18, < Af sp j >. This value is then used 
as the average of a poissonian distribution, that we assumed 
to assign a randomised number of spiral N sp i to each cluster, 



thus to determine the total cluster galaxies: A^oo = N'. 



R200 
gal 



spi • 



This formalism led us to derive an average fraction of about 
30% spiral galaxies in our sample of SDSS-maxBCG clusters 
(see TableQ]). The relation between N ga \ and the fraction of spiral 
galaxies inside the whole cluster is shown in Fig. El When con- 
sidering galaxies with the same magnitude limit (dashed line), 
the spiral fraction can be higher than 0.5 for groups and small 
clusters, and then lessens with increasing A^ai down to ^0.4 
for No a \=10. When considering all the galaxies included in our 
model (i.e. M, < — 16 for early-type andM r < — 18 for late-type, 
solid line) the spiral fraction becomes l ower by abou t 0.2. 

The morphology-density study of iBalogh et al.l (2004) be- 
comes useful also to assign a magnitude and a color to each spi- 
ral galaxy. For a given value of £5 it is possible to construct 
a probability distribution ¥(M r ,u — r) for a late-type galaxy to 
have a given magnitude M r and a given intrinsic color u — r. The 
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Fig. 2. Fraction of spiral galaxies within /?200 (see Sect. 12.1 
for the definition of /?2uo) as a function of Ng&. The solid line 
represents the result obtained when considering only galaxies 
wi th M r < —20 (i.e. the magnitude limit used to calculate £5 
in IBalogh et al.l 120041) . The dashed line is the result obtained 
by considering all the galaxies included in our model: late-type 
galaxies with M, < —18 and early-type galaxies withM,- < — 16 
from the SDSS-maxBCG catalogue. These functions have been 
calculated neglecting the evolution of the spiral fraction with 
redshift. 



former value is used to determine the luminosity L, in the r-band 

as 



-0.4M,. 



(7) 



where Lo,r=2.15x 10 34 erg/s is the luminosity of an object with 
M r = 0, while the latter is used to assign the spectral type, ac- 
cording to the intervals shown in Table [TJ To define these in- 
tervals we calculated the u — r colors of our GRASIL templates 
(2.40, 2.01 and 1.55 for the Sa, Sb and Sc, respectively) and 
defined the limits in order to associate to every galaxy the tem- 
plate that best approximates its value of u — r. We also checked 
the consistency of these num bers with the observ ed colors of lo- 
cal galaxies as reported by Fukugita et al .] (119951) for the SDSS 
photometric systerrfl. For what concerns starbur st galaxies, the 
limit o f u — r < 1.0 has been taken directly from Fukugit a et al.l 
(119951) for irregular galaxies, since the u — r of our M82 tem- 
plates (u — r — 2. 10) is not representative of the starburst galaxies 
population. However, we will show that the choice of this limit 
has a negligible impact on our final results. Finally, the luminos- 
ity Lir is obtained from Eq. [3] where Rirj- is chosen according 
to the spectral type. 

3.3. Evolution with redshift of spiral galaxies 

Despite being the SDSS-maxBCG a catalogue of relatively 
nearby groups and clusters (0.1 < z < 0.3), it is necessary to 
take into account the evolution with redshift of the late-type 
galaxi es properties with respect to the results of Balogh et al. 
(120041) which are derived from a local (i.e. z < 0.08) sample 
of SDSS galaxies. For the purpose of our modelisation we 
consider the sample of Balogh et al. (2004) as representative 



2 Referring to their Table 3, they obtain u — r colors of 2.26 for Sab 
and 1.68 for Sbc, thus close to our assumptions for the Sa/Sb and Sb/Sc 
limits, respectively. 
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of the galaxy population at zb — 0.04, the central value of the 
interval, and apply two independent evolution effects. 

Spiral fraction evolution 

It is generally accepted that in a hierarchical structure formation 
scenario spiral galaxies tend to evolve towards SO or ellipticals. 
This can happen either passively, with the consumption of the 
gas reservoir and its ejection with SNe explosions, or it can also 
be triggered by environ mental interactions with other galaxies 
and with the ICM (see iBoselli & Gava zzi 2006a, for a review): 
it is therefore expected that the fraction of late-type galaxies 
increases with redshift (the so ca lled Butcher-Oemler effect, 
iButcher & Oemlerl |1984l) . In fact, lLagana et al.1 <t.2009) find a 
change from about 0. 1 to 0.3 in the spiral fraction of 20 galaxy 
clusters in the interval < z < 0.25. Fitting the data on the 
fraction of star-forming galaxies of these clusters, they quantify 

its redshift evolution as: = 1.3 ±0.6 (pri vate communi- 
cation). A similar result has been obtained by iDe Lucia et al.l 

(120071) . who find ^ ~ 1.1 in the range 0.4 < z < 0.8. Given 
the redshift interval of the S DSS-maxBCG catal ogue, we take 
as a reference the value of lLagana et al.1 d2009h . Despite this 
sample contains rather massive objects (Icq T from ~ 3 to ~ 10 
keV), it should provide a first-order representation of the spiral 
fraction evolution in the redshift range of the SDSS-maxBCG 
catalogue. We therefore correct the number of spiral galaxies 
obtained in Section [3~2l with the following formula: 

<N spi (z) >= [1 + 1.3 (z-zb)} <N spi > z=0 ■ (8) 

Including this effect in our modelisation raises the global 
spiral fraction of our sample from 30% to 35%, as shown in 
Table Q] Alth ough this number may seem high for clusters (see 
e.g. lBai et al.l l2006. for the Coma cluster) we must consider that 
the SDSS-maxBCG catalogue is dominated by groups and small 
clusters, as mentioned in Section |2~T1 However, when applying 
our m odel to the case of MS 1054-03 as observed by iBai etall 
(120071) (z = 0.83, N 2 oo = 144), we obtain / spi = 0.16 which is 
consistent within a la limit with the value of 0.13 ±0.03 found 
by the authors. 

Luminosity evolution 

The current hierarchical structure formation scenario predicts 
that the peak of the star-formation is to be placed at z ~ 2 — 3 
and that at later epochs galaxies are mainly consuming their 
reservoir of gas. For this reason, when modeling the IR emission 
in the redshift range of the SDSS-maxBCG cluster sample we 
must consider that higher SFRs are expected with respect to the 
local galaxies. This effect has been widely observed and quanti- 
fied. We tak e as a reference t he ev olution of the IR-luminosity 
observed by iLe Floc'h et al.1 d2005l) : when parameterising it in 
the form of L iR (1 they obtain 0Cl = 3.2+Q2 over a 

sample of 2 635 objects identified at 24 /jm in the Chandra Deep 
Field South in the redshift range < z ^ 1 . This result has been 
obtained by analysing field galaxies and this might introduce an 
environment bias on our results. Anyway, since it is expected 
that local star-forming galaxies inside clusters were the latest 
accreted, we can assume that this evolution scenario is valid 
also in cluster and group envir onments. Moreover, as shown 
in figure 3 of iGiard et al.l (120081) . this assumption is backed-up 
by the Spitzer observation of the distant cluster MS 1054-03 
(IBai etal.ll2007h . Another bias may be introduced by the fact 
that these objects were identified directly in the mid-IR and 
not in the optical as for the galaxies of our sample. Therefore 
the galaxies of this sample are brighter than the ones of the 



SDSS-maxBCG, both in the IR and in the optical. Howev er, 
the total luminosity of the sample of ILe Floc'h et al. 1 (12001 is 
dominated by luminous (Lir> lO 11 !,©) massive galaxies, which 
are the ones that provide most of the signal also in the IRAS 
bands. Therefore, we assume that their result can be applied 
also to the SDSS-maxBCG cluster sample to describe its global 
IR luminosity evolution. Consequently, we determine the IR 
luminosity Lfr as a function of redshift: 

Lir(z)-Lir.o(1+z-zb) 3 - 2 , (9) 

where Lir is the luminosity obtained by the magnitude M r as 
from eq|7] 

We do not include any specific effect of spiral population 
evolution towards bluer types with increasing redshift. In fact, 
the fraction of the different spectral types of spirals in the red- 
shift range we are interested in is not expected to differ signif- 
icantly with respect to the loc al galaxies inside galaxy clusters 
(see e.g. IDe Lucia et al1 l2007). Moreover the consequent effect 
of slightly increased global luminosity is anyway partially in- 
cluded when assuming the luminosity evolution with redshift 
aforementioned. 

3.4. Uncertainties on the model parameters 

The choice of the values of the different parameters introduces 
a degree of uncertainty in our modelisation, the most important 
of which are connected with the choice of the redshift evolution 
parameters and the u — r color intervals used to define the spi- 
ral types (see Table [TJ. In order to quantify these uncertainties, 
we introduce two other sets of parameters to define a conser- 
vative and an extreme scenario and we will use the values of 
the fluxes estimated with these models to define the degree of 
confidence of our results. More precisely, for what concerns the 

redshift evolution we adopt the ±la v alues for -tff- = (0.7, 1 .9) 
and <Xl = (3.0,3.9), a s measured by Lag ana et alj d2009l) and 
ILe Floc'h et al.l d2005l) . respectively. Moreover, we also artifi- 
cially add a shift in the u — r color intervals of ±0.1 in or- 
der to obtain a redder (conservative scenario) and a bluer (ex- 
treme) spiral population. Anyway, these changes introduce only 
minor modifications in the galaxy population: the difference in 
the global spiral fraction is of ~ 2%, while the only significant 
change is in the fraction of starburst galaxies that goes from 
0.2% in the conservative model to 1.1% in the extreme one. On 
the whole, the difference in the final fluxes is mainly introduced 
by the modified luminosity evolution. The results on the galaxy 
population for the conservative model are also shown in TableQ] 

3.5. The IR luminosity function of the SDSS-maxBCG 
galaxies 

We show in Fig. [3] the luminosity function of the galaxies in- 
cluded in our reference model, for the different spectral types. 
The luminosity considered is given in the rest-frame 60 /jm 
IRAS band. 

The luminosities go from a minimum of - 10 8 L o for the 
fainter E/S0 galaxies up to a maximum of some lO n L for 
the brightest spirals. Each type of spiral has a luminosity range 
of about 2 orders of magnitudes: this is a consequence of our 
method (see Sect. 13.21) . that held the spiral magnitudes within 
the interval —23 < M r < — 18 which corresponds to a factor of 
1 q0.4 AM r _ 100 in luminosity with the only possible modifica- 
tion introduced by the redshift evolution. Anyway, we don't ex- 
pect this limitation to affect significantly our results. In fact, this 



M. Roncarelli et al.: IR properties of the SDSS-maxBCG clusters 



7 




8 9 10 11 

L°g(L [40 -60] „m/L s ) 

Fig. 3. Luminosity functions of each spectral type implemented 
in our model in the IRAS 60 /urn (i.e. top hat band-pass within 
[40,80] fjm). Spectral types are color-coded as follow: (i) early- 
type galaxies as (dark-red) for the BCGs and (light-red) for the 
other E/SOs (we assigned to every galaxy the average luminosity 
of the E/SO members of each cluster, as shown in eq. |4j; (ii) nor- 
mal late-type galaxies, with Sa, Sb, Sc going from (light-blue) 
to (dark-blue), respectively; (iii) starburst galaxies as (green). 

limit range is wider than the one covered by the distribution of 
the BCGs, which are the only galaxies for which we have di- 
rect observational constraints. On the contrary the E/SO curve 
appears narrower because we considered the average luminosity 
of the E/SO galaxies inside each cluster, thus limiting the disper- 
sion of the E/SO galaxies to the dispersion of the cluster lumi- 
nosities. 

When considering the luminosity function in the 100 fjm 
IRAS band the range and shape is similar for all the spectral 
types. Only starbursts are a factor ~3 fainter in this band due to 
the drop in luminosity (see Fig. [TJ which is due to an average 
higher galactic dust temperature compared to non-starburst. 

4. Other galactic components 

In this Section we try to put some constraints on the impact 
of three components of galactic origin which are not included 
in our modelisation of the IR emission of cluster galaxies: a 
missed population of faint galaxies, the possible presence of 
dust-embedded AGNs and the IR emission coming from heavily 
obscured star-forming galaxies. In all of these cases we conclude 
that their impact on our final results is probably very small, if not 
completely negligible. 

4.1. Faint galaxy population 

As mentioned in Sect. [3] our modelisation includes red galax- 
ies with M r < -16 (the limit of the SDSS-maxBCG catalogue) 
and blue galaxies with M r < — 18, thus neglecting the signal of 
fainter objects, that may not be identified by the SDSS observa- 
tions although present inside the cluster. We try to obtain a rough 
estimation of the impact of these objects on the tota l emission by 
referr ing to the luminosity functions obtained by iBaldrv et al.l 
(12004 (see their Fig. 7). 

For what concerns early-type galaxies, the faint-end of the 
LF (M* = -21.49, a = -0.83 in the Schecter function param- 



eterisation) has a negative slope, so the number of galaxies is 
expected to diminish at higher magnitudes. For this reason no 
significant impact can be associated to these objects. 

On the other side, the faint-end of the blue galaxies LF 
(M* = —20.60, a = —1.35) has a positive slope, so an increas- 
ing number of galaxies is expected with lower luminosities. For 
simplicity, we consider the assumption that all galaxies with 
M r < — 1 8 have been included in our modelisat i on, th en we in- 
tegrate the late-type LF §b(M r ) of IBaldrv et alJ d2004 splitting 
it in two at this magnitude limit. We obtain 

^faint _ ri%^{M r )L{M,-)dM _ q ^ 

^bright Jl23<S?b(M r )L(M r )dM ' ' 

where L(M r ) is the luminosity as a function of the magnitude M r , 
thus indicating that the emission of the faint galaxy population 
is marginal with respect to the bright one. 

Moreover, since the LF of lBaldrvetall d2004l) is obtained 
with the observation of field galaxies, it is reasonable to expect 
that in cluster environments the faint-end of the LF would be 
shallower, if not even negative, due to the processes of merging 
that affect particularly the smaller objects. Therefore the ratio of 
0.13 obtained in eq.[l0]can be safely considered an upper limit 
of the true value. For these reasons, and considering the assump- 
tion of the connection between the luminosity in the optical and 
IR bands used throughout this work, we can conclude that even 
if we can not exclude the presence of the signal of a popula- 
tion of faint (M r > — 18) unresolved galaxies, we expect that the 
contribution of faint star-forming galaxies to be marginal. 

4.2. Dusty AGNs 

When estimating the total IR luminosities of the galaxies of our 
sample, we do not take into account what could be the contribu- 
tion of AGN deeply embedded in dusty cocoons within cluster 
galaxies. 

In fact, recent results indicate the existenc e of a population 
of heavily absorbed A GNs in the field (see e.g. lFiore et all 2009: 
lLanzuisi e t al. 2009]): these objects are detectable at IR wave- 
lengths. Although optical observations indicate a very small 
fraction of cluster g alaxies with detected AGNs (~ 1%, see e.g. 
iDressler et al.ll 19991) . some X-ray observations have detected an 
excess of point sources associ ated to AGNs in cluster fields. For 
example, Mart ini et alJ J2006) observed with Chandra a sample 
of 8 low redshift (z < 0.3) clusters, finding ~5% of their galaxy 
members hosting an AGN, most of which are not detected with 
optical surveys. These AGNs are also present in E/SO galaxies, 
and whilst they are obscured in the optical (their emission being 
heavily absorbed), they could be bright in the IR. Therefore, they 
could contribute to the total IR cluster emission. 

It is difficult to quantify with precision the impact on the to- 
tal IR signal of t hese kinds of obj ects. We take as a reference 
the work done by iBai et al . (2007), who studied the IR proper- 
ties of the galaxies of MS 1054-03 (z = 0.83) with the Spitzer 
satellite. They identified eight point sources with X-ray and ra- 
dio observations that could be associated with AGNs: anyway, 
only three of these objects have an IR counterpart over the 144 
IR-detected cluster members, and only one of these is associ- 
ated with a bright star-forming object. For these reasons the au- 
thors conclude that the contamination f rom AGN s is neg ligible. 
Similar conclusions have been drawn by Ba i et al. (2006) in then- 
study of the Coma cluster (z — 0.023). 

Even if they are based on single cluster studies, these results, 
which have been obtained both at redshift higher and lower than 
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our cluster sample, indicat e that the impact of AGN contamina- 
tion in the observations of iGiard et al.l (120081) is probably very 
small. 



4.3. Heavily obscured star-forming galaxies 

Optical observations have highlighted the fact that some galaxy 
clusters host heavily obscured star-forming galaxies. In the 
framework of our modelisation, the presence of these objects 
would lead to an u nderestimate of thei r IR emission based on 
the optical one (see lMetcalfe et al1l2005l for a revi ew). In partic- 
ular, IR observat ions on A1689 (iDuc et al.ll2002l). J 1888.16CL 
dDuc et al.ll2004h . C L0024+1654 dCoia et al.l l2005h and, more 
recently, on A1758 dHaines et al.l 120091) have revealed galax- 
ies with much higher SFRs with respect to what is expected 
from optical diagnostics (e.g. [Oil]). However, all of these clus- 
ters show clear signs of recent dynamical activities, like ma- 
jor mergers and significant accretion of galaxies from the field, 
thus suggesting that these phenomena are responsible of the star- 
formation triggering. Since our cluster sample is at a relatively 
low redshift, we do not expect it to contain a high fraction of 
dynamically active haloes and, therefore, to be significantly af- 
fected by the presence of heavily obscured star-formation. 



5. Reconstructing the stacked IR flux 

As mentioned before the main objective of this work is to ex- 
plain the origin of th e IR emission observed in t he direction 
of gal axy clusters by Montier & Giard (2005|) and IGiard et al. 
(2008). Therefore, we want to compare the flux expected out of 
our model from the galactic emission in the 60 ^m and 100 /jm 
IRAS bands with the one measured by Giar d et all J2008). In or- 
der to do so, we need to compute for every single cluster the 
expected flux taking into account the instrumental beam and the 
spectral band pass of the IRAS satellite. 

We assume that every galaxy is placed at the redshift z of 
the cluster to which it belongs and we compute its flux in a 
given band 



possible loss of signal from galaxies distant from the cluster cen- 
ter due to the convolution with the instrumental beam. In order to 
model this effect, we have to distribute the cluster galaxies in the 
cluster potential well. We thus randomly assign to every galaxy 
a distance from the BCG (considered to be at the cluster center, 
thus at the center of the 10' field) by assuming that their spa- 
tial distribution follows a NFW profile dNavarro et al TT997 ). The 
two parameters needed to characterise the NFW profile are R200 
and the concentration c. They were obtained from the richness 
^gai 00 °f eac h ( see Sect. 12. It and by adopting the N - M 
scaling relation of lRvkoff etaTI (120081) derived from the SDSS- 
maxBCG, and the (M,z) - c relation of iDolag etakl d2004l) de- 
rived from numerical simulations. Then for every galaxy we 
compute the fraction of the signal that falls inside the observed 
region in the two bands. This effect proves to be completely neg- 
ligible as globally less than one per cent of the signal falls out- 
side the field in both bands. Since galaxies of later spectral types 
are expected to be more spread than ellipticals, as a result of 
their recent accretio n into the cluster (see, e.g., the discussion in 
Pope sso et al.l 12005), we repeated our estimation on the possi- 
ble lost signal by assuming that late-type galaxies are distributed 
uniformly in a sphere of radius /?20o: even with this extreme hy- 
pothesis the amount of signal that is expected to fall outside the 
field remains negligible (about 3% of the total). The low impact 
of this effect on the global results can also be seen from the dis- 
tributions of the angular sizes of /?200 obtained with this method 
shown in Fig.|4]for different redshift bins. More than 70% of the 
clusters have angular sizes lower than 5' and only in the lowest 

redshift bin some objects (43) exceed the 1 0' aperture. 

With the stacking technique adopted bv lGiard et al.1 (12008). it 
is clear that also the signal coming from foreground/background 
objects is present in the observed fields. Anyway, as said in 
Section [2] the authors adopted a background subtraction based 
on the signal in the fields adjacent to every map, relying on the 
high statistical robustness of their sample. For this reason we can 
safely conclude that the emission of foreground/background ob- 
jects has been successfully subtracted, thus we do not need to 
include it in our modelisation. 
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is the luminosity distance of the cluster (Hq = 100/i km s 1 
Mpc -1 ), (A,) is the IRAS spectral response function in a given 
band (i.e. A-o = 60, 100 /jm) which is convolvecQ with the SED 
L(X) of the spectral type of each cluster galaxy, normalised ac- 
cording to the value of Lir (see Sect. 

The stacked signals by IGiard et all (120081) are integrated 
fluxes over an angular area of 10' of radius centered on the clus- 
ter (see Sect. [2] for details). Since the IRAS FWHM beams are 
4' and 4.5' at 60 and 100 /jm respectively, thus comparable to 
the size of the observed field, we need to take into account the 



This quantity corresponds to a flux per unit frequency (e.g. Jy). 
Anyway we prefer to use the notation F^, rather than F v , to refer directly 
to the fluxes F^ and Fjoo in the 60 and 100 /jm bands, respectively. 

4 For these calculations we used the too ls of the DUSTEM c ode, an 
updated version of the emission model of lDesertetai] d!990h . which 
contains details on the IRAS response functions in the different bands. 



6. 1. Predicted and observed IR fluxes 

Our estimations of the 60 and 100 /jm fluxes from cluster galax- 
ies are reported in Table [2] associated to the emission calculated 
with our model for the different galaxy populations. 

The total flux due to the early-type galaxies estimated with 
our model is 39.3 Jy at 60 /jm band and 79.4 Jy in the 100 
band, accounting for about 7% of a nd 6%, respectivel y, of the 
fluxes measured in the same bands bv lGiard et al.l(l2008l) . In both 
bands about 20% of the E/S0 signal comes from the BCGs. 

As expected, late-type galaxies constitute by far the most sig- 
nificant contribution to the IR emission. According to our refer- 
ence model they contribute to 645.2 Jy at 60 and 1904.8 Jy 
at 100 /jm thus accounting for about 95% of the total galactic 
emission. This contribution comes mostly from the Sb and Sc 
population, with the Sc being widely dominant in the 100 /jm 
band. 

Despite their high luminosities, we predict that starburst 
galaxies do not provide a significant contribution to the IR emis- 
sion, due to their low expected number. We obtain a contribution 
of 15.4 Jy at 60 and 12. 1 Jy at 100 /jm corresponding respec- 
tively to about 2% and 0.5% of total predicted signal. Only in 
our extreme scenario their contribution becomes non-negligible 
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Fig. 4. Distribution of the angular sizes of /?200 for the cluster sample in 4 redshift intervals chosen in order to contain the same 
number of clusters (i.e. 1869). The bin size of each histogram is 0.5 arcmin. 



Table 2. Estimated fluxes in the 60 /jm and 100 /jm IRAS bands 
from the galaxy population of the SDSS-maxBCG clusters. 



Spectral type 




F\00 




(Jy) 


(Jy) 


-EaHy-tvpe 






BCG 


8.2 


16.5 


E/S0 


31.1 


62.9 


(Total early-type) 


(39.3) 


(79.4) 


-Late-type 






Sa 


36.4 


58.6 


Sb 


259.5 


621.3 


Sc 


333.9 


1133.4 


Starburst 


15.4 


12.1 


(Total late-type) 


(645.2) 


(1825.4) 


Total 


684.5™ 


1904.8+^'j 


Observed 


570.1±36.1 


1359.9±249.1 



We show the total fluxes and contribution of each galaxy population 
included in our model. The quoted errors on the total values correspond 
to the differences between our reference model and the conservative 
and extreme scen arios described in Section [3^41 The measurements by 
iGiard et al] i2008h are reported together with their la error bars. 



in the 60 /jm band, reaching the 5% of the total predicted flux. 
This low contribution agre es with the low rate of s tarburst galax- 
ies as found in the field bv lLe Floc'h et alJ (2005) at the redshift 
range of the SDSS-maxBCG catalogue (i.e. 0.1 < z < 0.3). 

The total fluxes associated to the galactic emission predicted 
by our reference model are 684.5 [546.5,895.8] Jy at 60 /jm 
and 1904.8 [1475.7,2521.9] Jy at 100 ^m (the bracketed interval 
indicate the values derived from our conservative and extreme 
models, see Sect. l3.41 >. It appears that the reconstructed IR emis- 
sion due to the galactic dust emission of t he cluster member s 
can explain the entire signal measured by IGiard et alJ J2008), 
with an indication that our reference model overestimates the to- 
tal flux, particularly at 100 /jm. Indeed, these authors obtained 
570.1 ±36.1 Jy and 1359.9 ± 249. 1 Jy at 60<um and 100 <um, re- 
spectively. We will propose an explanation of this discrepancy in 
the next sections. However, when considering our conservative 
scenario, the predicted emission is in good agreement with the 
total measured signal in both bands. 

Given these results, modulo the uncertainties of our modeli- 
sation, we obtain that the IR emission of the galaxy members is 
consistent with the total observed emission of our clusters sam- 
ple, leaving little space to the possible presence of other compo- 
nents like intracluster dust. 
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Fig. 5. Average cluster luminosities (see eq. [T3J as a function 
of redshift. Diamond s represent the me asurements and associ- 
ated la error bars by iGiard et al.l (I2008I) in 4 redshift bins, de- 
fined in order to contains the same number of clusters (i.e. 1869). 
The solid line represents the average luminosity of the clusters 
in our reference model. The two dashed lines show the corre- 
sponding luminosity of the conservative and extreme model (see 
Section [3~4l for the details). 



6.2. Redshift evolution 

In this Section we characterise the redshift evolution o f the av- 
erage clusters luminosity. Following D esert etafl (Il990h . we de- 
fine the total luminosity in the 60 and 100 fjm IRAS bands as 



100 : 



- 4nd? x 



Av Av 
X—F x {60fim) + X—F x (l00fmx) 



(13) 



where Av and AX, are the bandwidths in frequency and wave- 
length, respectively, of the two IRAS bands. This formula repre- 
sents a good approximation of the total luminosity of our haloes 
at X > 40^m. 

We show in Fig.|5]the average luminosity as obtained from 
eq.Q~3]for the clusters at differe nt redshifts, both as predicted by 
our model and as observed by IGiard et all (120081) . Our 4 bins 
are defined in order to contain the same number of clusters (i.e. 
1869). This analysis shows clearly that the discrepancy between 
our reference model with respect to the observed fluxes is due 
to the low-redshift clusters only, which are responsible for most 
of the predicted flux. In fact, while our reference mode l is com- 
patible within la with the results of IGiard et "all d2008l) in the 3 
high-redshift bins, in the first bin (0. 10 < z < 0. 17) the predicted 
luminosity is higher at 6a confidence. Our conservative scenario 
is also compatible with observed luminosities at z > 0.17, while 
the discrepancy for the low-redshift clusters persists at more than 
3 a. On the contrary, our extreme model is excluded by the results 
of lGiard et all d2008l) through all the redshift range. 

6.3. Selection bias 

Our results indicate that although our model well describes the 
global IR emission of the SDSS-maxBCG galaxy clusters at 
z > 0.17, the predicted galactic emission at low redshift clearly 
exceeds the measured stacked signal. The strong decrease in the 
IR luminosity is present also in the N ED and NSC clusters sam- 
ples analysed by IGiard et all d2008l) and, as discussed in their 



Fig. 6. Average richness as a function of redshift. Diamonds rep- 
resent the average values of A^ 00 (i.e. the number of E/S0 mem- 
bers) with associated la errorbars in 4 redshift bins defined as 
in Fig.[4]and[5] Triangles represent the values of A200 (i-e. total 
cluster members, see the definition in Section [3~2T i. 



work, it can be interpreted as a selection effect which biases to- 
wards rich and massive clusters at high redshift as it is evident 
for the NSC sample (see their Fig. 6). However, this is not the 
case for the clusters object of our analysis. In Fig.|6]we show the 
average cluster richness in the same redshift bins both for the 
values of A™ , directly obtained from the SDSS-maxBCG cat- 
alogue, and A^oo which includes late-type galaxies, as described 
in Section 13.21 It is clear that no selection bias is present in 
our clusters sample: on the contrary low-redshift clusters show 
slightly higher values of A™ . Even when including late-type 
galaxies this trend does not change significantly. 

6.4. Effect of the cluster environment 

These considerations indicate that the discrepancy between our 
model and the decrease of luminosity towards lower redshift 
must be connected with an evolution of the galaxy IR luminos- 
ity driven by the cluster environment, associated to gas and/or 
dust removal and consequent star-formation quenching. This 
can happen via ram-pressure stripping and tidal interaction. 
This picture is in agreement with the low number of high IR- 
to-optical galaxies obser ved in local galaxy clusters (see, e.g., 
Bicav & Giova nelrilll987l) . A similar phenomenon is also seen 
in the obse rved properties of th e galaxies of the Virgo clus- 
ter. In fact, iBoselli et al.l (2006) observed a truncation in the 
disk of NGC 4569 both in the UV and in the IR (8-70 /urn) 
that they associate to a r am-pressure st r ipping of the external 
regions. More in general, iGavazzi et al.l (|2006) observe a con- 
sistent H a deficiency in Virgo galaxies which indicates a low 
SFR a nd diminished dust heating /emission (see also the discus- 
sion in lBoselli & Gavazzil2006al) . In this scenario, once the star- 
formation quenching happens, the SEDs of the galaxy members 
would be consistently modified at the IR and UV wavelengths, 
being not compatible anymore with field galaxies, thus explain- 
ing the excess predicted by our model. However, we must point 
out that the environmental effects on the Virgo cluster are likely 
to be much stronger than what is expected for the population 
of clusters considered in our analysis. In fact, the Virgo clus- 
ter is about one order of magnitude more massive than the ma- 
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jority of the SDSS-maxBCG galaxy clusters and it shows sev- 
eral evidences of recent dynamical activity. On the contrary, be- 
ing the SDSS-maxBCG cluster sample mainly constituted by 
small haloes (see Section 12. U . in the context of the hierarchi- 
cal structure formation scenario it is reasonable to expect that 
they correspond to relaxed and dynamically old systems (see, 
for instance, the results on the con centration-mass relation of 
iDolae et al.ll2004tlBuote et alJl200l . 

6.5. Constraints on intracluster dust 

If we take as a reference our conservative scenario and com- 
pare its results with the average observed luminosity in the three 
high-redshift bins, we can estimate an upper limit on the possi- 
ble emission due to extragalactic dust of about 10% of the to- 
tal luminosity, which translates into a dust- to-gas mass abun - 
dance of Z<] < 5 x 10~ 5 (see the discussion in lGiard et al.l l2008). 
This figure is in agreement with the expectations from the- 
oretical models (Popescu et aT] l2000l lAguirre et al]l200ll) and 
with current observatio nal upper limits on the dust abundance. 
IChelouche et al.l d2007l) analysed the reddening of a sample of 
quasars in the direction of the clusters of the SDSS-maxBCG 
catalogue, obtaining an es timate of Za g 10~ 5 . A similar result 
has also been obtained by iBovv et ail (120081) by measuring the 
dust absorption on the spectra of galaxies located behind local 
(z ~ 0.05) galaxy clusters. 

If on one side we do not see an evidence of intracluster dust 
emission in our sample, these last considerations leave an open 
question of how much dust has been lost by cluster galaxies pol- 
luting the ICM and if it can live enough to produce a signif- 
icant diffuse IR emission. According to Popescu et al. (2000), 
the dust stripped from infalling galaxies will probably remain lo- 
calised close to their parent galaxies without diffusing efficiently 
into the ICM. In any case, the strong d e crease in the total clus- 
ter luminosity observed by iGiard et al.l (120081) compared to our 
results, suggests that, if present, the signal of intracluster dust 
should be very small. 

7. Summary and conclusions 

In this work we performed, for the first time, a thorough modeli- 
sation of the overall IR emission of galaxy clusters due to cluster 
galaxies IR emission. We tested the results of our model against 
the statistical stacking measurements of clusters IR emis sion in 
the 60 and 100 IRAS bands bv IGiard et all d2008l) . mak- 
ing use of the SDSS-maxBCG catalogue of groups and clusters 
dKoester et al.ll2007l) . 

We used the available SDSS-maxBCG data on the lumi- 
nosity in the r-band of the cluster early-type members, and 
converted them into IR luminos ities by adopting the model 
templates of the grasil code dSilvaet all 1 19981) . Since the 
SDSS-maxBCG catalogue does not contain information on 
the late-type galaxies component, we used the resu l ts on the 
morphology-density relation obtained bv lBalogh et al.ld2004l) on 
a local SDSS galaxy sample to construct a model to statistically 
associate to every cluster its spiral galaxy population, namely 
the number of spiral members, their spectral type and their M r . 
Again, this information on the optical properties has been used to 
obtain the corresponding IR emission by making use of 4 other 
GRASIL templates to represent normal late-type (Sa, Sb and Sc) 
and starburst galaxies. We also included in our model the ex- 
pected redshift evolution of the spiral fraction and the IR lumi- 
nosity by using results on observed galaxies in order to account 
for the possible difference in the spiral galaxy population of the 



SDSS-maxBC G catalogue (0.1 < z < 0.3) and the more local 
objects of the lBalogh etai] d2004 sample (z < 0.08). 

Finally, we used our predictions to calculate the total ex- 
pected galactic flux, by considering also the possible loss of sig- 
nal due to the IRAS beam smoothing and the IRA S response 
functio n, and compared it with the measurements of Giard et ail 
d2008l) . 

Our main results can be summarized as follows. 

i) According to our model, late-type galaxies represent about 
~ 35% of the galaxy population in nearby groups and clus- 
ters (i.e. z < 0.3). This figure is coherent with the fact 
that this sample is dominated by groups and small clus- 
ters (NgJ 00 = 10 - 20), thus letting this fraction of late-types 
falling between the field (/ sp ; 0.5) and the rich cluster 
(/s P i « 0.15) regime. 

ii) As expected, normal late-type galaxies constitute the most 
important contribution to the total IR luminosity at 60 and 
100 /jm accounting for ~ 95% of the total galactic emission. 
Oppositely, the impact of starburst galaxies is marginal. 

iii) Early-type galaxies dominate the faint end of the LF in the 
IR and they account for the remaining ~ 5% of the galactic 
contribution of the cluster IR emission. 

iv) Our model shows that the total flux estimated from the 
galaxy accounts, wi thin uncerta i nties, for the entire stacked 
signal measured bv IGiard et al.l d2008). With our reference 
model we obtain 684.5 Jy at 60 ^um and 1904.8 J y at 100/im 
which exceeds the stacked fluxes measured by IGiard et ail 
d2008l) . However, when considering the uncertainties in our 
parameters, we showed that in a conservative scenario this 
discrepancy disappears. 

v) We compared the redshift distribution of the I R emission 
predic ted by our model to the measurements of IGiard et al.l 
(2008) and found that the excess flux is present only in the 
clusters at lower redshift (z < 0.17). 

The results presented here show that the IR emission of 
galaxy clusters can be explained with its galactic component 
only, leaving very little room to the possible presence of any dif- 
fuse emission associated to intracluster dust. In this framework, 
the upper limit on the dust-to-gas mass abundance obtained by 
IGiard et alJ d2008) can be reduced by an order of magnitude, 
down to Zd < 5 x 10~ 5 . This result is in agreement with current 
estimations on the dust abundance in the ICM obtained from ex- 
tinction and reddening measurements in the dire ction of SDSS 
clusters dChelouche et al.ll2007l: IBovv et alJl2008l) . 

The lack of diffuse dust, however, does not mean that 
the environment is not influencing the IR properties of its 
galaxy members. The fact that the predicted excess emission 
is concentrated in the lowest redshift objects indicates that 
the IR emission of late-type galaxies in local clusters is not 
completely compatible with their field equivalents. This result 
is in agreement with the lack of bright IR galax i es obs erved 
in several local clusters (see iBoselli & Gavazzil [2006a, and 
references therein) and it is likely connected to the quenching 
of the star-forming activity driven by the cluster environment 
on its galaxy members, due to gas or dust removal: in this last 
case, the dust injected in the ICM is probably quickly depleted 
via sputtering processes. However, it is not clear why this could 
act significantly only at z < 0.2 when most of the large-scale 
structures are already formed. For this reason, and given the 
limited redshift range of the SDSS-maxBCG catalogue, it would 
be interesting to extend our analysis to higher redshifts, where 
the SFR is expected to be higher and more dynamical activity is 
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